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ABSTRACT 

The origin of Earth oceans is controversial. Earth could have acquired its water either from 
hydrated silicates (wet Earth scenario) or from comets (dry Earth scenario). [HD0]/[H20] 
ratios are used to discriminate between the scenarios. High [HD0]/[H20] ratios are found 
in Earth oceans. These high ratios are often attributed to the release of deuterium enriched 
cometary water ice, which was formed at low gas and dust temperatures. Observations do 
not show high [HD0]/[H20] in interstellar ices. We investigate the possible formation of high 
[HD0]/[H20] ratios in dense (nn >10^ cm-^) and warm gas (T =100-1000 K) by gas-phase 
photochemistry in the absence of grain surface chemistry. We derive analytical solutions, tak- 
ing into account the major neutral-neutral reactions for gases at T > 100 K. The chemical 
network is dominated by photodissociation and neutral-neutral reactions. Despite the high gas 
temperature, deuterium fractionation occurs because of the difference in activation energy be- 
tween deuteration enrichment and the back reactions. The analytical solutions were confirmed 
by the time-dependent chemical results in a lO^'^ M0 disc around a typical TTauri star us- 
ing the photochemical code ProDiMo . The ProDiMo code includes frequency-dependent 
2D dust-continuum radiative transfer, detailed non-LTE gas heating and cooling, and hydro- 
static calculation of the disc structure. Both analytical and time-dependent models predict high 
[HD0]/[H20] ratios in the terrestrial planet forming region (< 3 AU) of circumstellar discs. 
Therefore the [HD0]/[H20] ratio may not be an unique criterion to discriminate between the 
different origins of water on Earth. 

Key words: Astrochemistry 



1 INTRODUCTION 



The water deuteration abundance ratio [HDO]/[H20] is often used to determine the temperature at which water was synthesised. Current 
sub-millimeter observations already provide measurements of [HDO]/[H20] for many astronomical objects includin g comets, hot-cores, an d 
protoplanetary discs. In particular, a high [HD01/[ H201 has been observe d toward the protoplanetary disc DM Tau jCeccarelli et al 
although the HDD detection remains controversial ( Guilloteau et al.lEoO( j). The Herschel Space Telescope will be capable to detect cool and 
warm water in discs while the Atacama Large Millimeter Array (ALMA) will spatially resolve HDD emission. 

The water deuteration abundance ratio [HDO]/[H20] in hot cores around massive young stellar objects is ~ 3 x lO"'^ 1 Gensheimer et al] 

1 199 ^. around a factor ~ 10 enhancement compared to the cosmic [D]/[H] value set at the Big Bang (~ 1.5x 10~^. lLinskvl2003 h. The enrich- 
ment suggests that the chemistry oc curs in a low density, col d medium (T < 100 K), where HDD synthesis begins with the deuteration of the 
molecular ion H+ into H2D+ (e.g.. lRoberts & Millaj|200oh . Further gas-phase reactions lead to a high [HD0]/[H20]. But at temperatures 
greater than 100 K, which is the typical gas temperature of hot cores, the initial deuteration reaction is inefficient and the [HDO]/[H20] is 
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close to the cosmic value. Therefore current chemical models need to invoke the evaporation of deuterium enriched water ice when the tem- 



perature reaches 100 K to explain the deuterium enrichment in water. The observed (HD0/H20)ice ratio is however too low dPartois et al 
2OO3I : IParise et allboOSh . A second difficulty arises from the high observed water ice abundanc e of 10 ^ -10 *, which is two orders of 



magnitude larger than the gas phase abundance of water in the hot core around IRAS 16293-2422 ( jParise et al.ll2005h . 

After decades of research, the origin of water on Earth remains a major subject of debate fe.g. jNutljbooi for a review). Two majors 
models exist. In the first model, the Earth accretes dry: the building blocks of the Earth are made of dry rocks, composed only of silicates 
and carbonaceous materials. Most of the water is brought afterward by comets during the phase of heavy bombardment. This model is 
supported by the similarity between the Earth Mean Ocean Water (SMOW ~ 1.49x 10"'') [HDO]/[H20] and the cometary value, although 
the later value seems too high. Another support to the "dry" model is the low [HD0]/[H20] predicted at temperature T > 100 K by 
thermochemical equilibrium models. In the second model, most of the water on Earth comes from the release of water vapor trapped inside 
the water-entrapped planetesimals like carbonaceous chondrites upon im pact or during volcanism. Water-rich planetesimals located at 2- 
3 AU are perturbed by the giant planets and collide with the young Earth jMorbidelli et al ] |200d l iGomes et al I booi l [Raymond et al ] |2004 



2005D . Those plan etesimals contain w ater in form of hydrated silicates. The average D/H in carbonaceous chrondrite is similar to the value 



in the Earth ocean jRobert et allboOOh . although the composi tion of carbo naceous chrondrite does not match the composition of Earth' crust 
I Righter et al. The Earth is said to have accreted "wet" iDrakel2005l) . iydrous material does not necessary need to be brought from the 

asteroid belt area (2-3 AU). Recent study on adsorption of water molecules onto fractal dust grains shows that hydrous material could have 
been present in the vicinity of 1 AU jstimpfl et allbood) . In this scenario water in the gaseous solar nebula sticks onto the silicate surface of 
fractal grains, which coagulate and grow into the Earth. The adsorbed water would reflect the deuterium enrichment of the gas-phase water at 
I AU. The major difference between the two "wet" scenarios is that water was present in the early-Earth in the latter scenario . Geochemical 
findin gs suggest that hydrosphere and continental crust were already present in the first ~ 600 million years of Earth history jHopkins et al 



2008h . The presence of water and also of organic matter at the formation of the Earth has profound astrobiology implications for the origin 



of life. 

One major weakness of the "wet" scenario is that the deuterium enrichment predicted by equilibrium chemistry is too low at temperatures 
greater than 100 K. However there are chemical routes to obtain high [HD0]/[H20] ratios at temperatures greater than 100 K in non- 
equilibrium chemistry. Photodissociation o f H2O and HDO followed by reformation at the surface of protoplanetary d i scs an d turbulent 
mixing can also change the [HDO]/[H20] 1 Yung et al. I988h. Vertical turbulent mix ing was invoked bv lUvons & Young (2005) to explain 
the oxygen anomalies in the early solar nebula. Alter natively Genda & Ikoma (2008) propose that the [HD0]/[H20] ratio on Earth changed 
during the evolution of the Earth from its formation. IWillacv & Woods C2009.) studied the deuterium chemistry of the inner disc around a 
TTauri star and found relatively high HDO abundance in the warm molecular layer. 

Four direct physical processes lead to isotopic fractionation in the gas phase. First, heavier isotopologues diffuse slower than the lighter 
main isotopologue. The diffusion results in zones of varying [HD0]/[H20] ratios. This process is negligible in turbulent media when the 
mixing timescale is much lower than the diffusion timescales. Second, due to their larger reduced mass, isotopologues usually have a lower 
zero point energy. The vibrational energy levels are located lower, which increases the density sum of energy and hence reduces the vapor 
pressure. Hence, the water ice/gas ratios in the interstellar medium is not at equilibrium. Third, heavier isotopologues are chemically more 
stable because the binding energy between O and D is stronger than between O and H in the water molecule. Fractionation reactions 
dominate at low temperature because the inverse reactions have a minimum energy barrier equal to the difference in binding energy. Finally, 
isotopologues differ in their photodissociation cross section due to changes in the selection rules. A combination of the last two processes 
together with chemical diffusion and turbulent mixing probably determine the deuterium fractionation in the inner region of protoplanetary 
discs where terrestrial planets form. 

Interstellar chemical models often neglect or include only few neutral-neutral reactions because their rates are small at temperatures 
below ~ 100 K but they become competitively fast with the ion-neutral reactions at a few hu ndred degrees. F or example, in dense and warm 
gases, the central species is the hydroxyl radical OH, which can react with H2 to form H2O (Thi & BM 20051) . 

In this paper, we explore neutral-neutral and photochemical reactions relevant to the formation and destruction of H2, HD, OH, OH, 
H2O and HDO at gas temperature between 100 and 1000 K. We first focus on steady-state abundances and perform an analytical analysis. 
Analytical analysis makes it possible to determine the main formation and destruction paths. After establishing the possible reaction paths, 
we use the photochemical code ProDiMo to compute a time-dependent chemical structure for a typical TTauri disk to check the validity of 
our assumptions in the analytical analysis. In the rest of the paper, the chemical reaction rates are given in Table|4l[5]|6l|7l and[8] In particular, 
the most important reactions are summarised in Tables ID 



2 STEADY-STATE ANALYTICAL SOLUTION 
2.1 Formation and destruction of H2O and HDO 



We consider a gas mixture at a single temperature, density, irradiated by UV photons. In thermochemical equilibr ium the deuterium exchange 
equilibrium reaction between water and molecular hydrogen (the most abundant reservoir of deuterium) is (e.g.. lRobert et aklboOOh 



HD + H2O ^ HDO + H2. 



(1) 
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The equilibrium constant of this reaction for gas between 273 and 1000 K is i Richet et all 19771) 



^ _ PhdqPh, ^ 0.22 X 10"^ 

" PhdPh.o - ^ ^ ' 

The pressure dependence of K is negligible. At equilibrium, the deuterium enrichment decreases quadratically with temperature. A gas 
mixture is in thermochemical equilibrium when each chemical reaction is exactly balanced by its reverse reaction. Dissociations upon 
absorption of UV photons are not taken into account. 

A gas is in chemical steady-state (also called stationary state) if the various abundances do not vary with time but the formation and 
destruction reactions for a given species are not necessarily the same. An analytical solution to the chemical network for the synthesis of 
H2O and HDO can be obtained by identifying the main formation and destruction chemical reactions. At high temperature (T >100 K) and 
dens ity, it is well established that water is p rimarily formed via the reactions between the hydroxyl radical OH and molecular hydrogen H2 
fe.g. lThi&Bikl2005 : iGlassgold et al horn . High temperature gases are needed to overcome the large energy barrier of the reaction (reaction 
5 with rate fcs, Ea=l660 K). The main destruction mechanism is likely photodissociation in low Ay regions. Other water destruction mecha- 
nisms are reactions with atomic hydrogen, and with ionised helium in shielded regions. We assume in this study that the dust temperature Td 
has equilibrated with the gas temperature. At Td > 100 K most water molecules remain in the gas phase. The rate of water formation is thus: 

^^^^ = (^5(^2] + kis[HD] + k7[OH])[OH] - (Js + kslH] + kig[D] + k4i[He+] + kcp,3)[H20], (3) 

where fcs, fcig, and k^ are the rates of formation of water via the reactions OH-1-H2, OH-l-HD, and OH-l-OH respectively, Jg is the photodisso- 
ciation rate, fcg is the water destruction rate by H, fccp,3 is the cosmic-ray induced photodissociation of water, and k^i is the rate of destruction 
by He"*". At steady-state (d[H20]/dt — 0), we obtain the ratio 

[H2O] ^ k4H2\ 

[OH] Je + fcs [H] + fc4i lHe+] + fc^p.s ' 

which will be used later in this paper. We have neglected the formation of water via reaction of OH with HD and OH and destruction of 
water via atomic deuterium because the rates are orders of magnitude small er than the rates of other paths. Deuterated water can be formed 
via the reaction of HD with OH as well as via OD -1- H2 ( Bergin et al.ll 19981) . Destruction occurs via photodissociation, reaction with atomic 
hydrogen and ionised helium. 

^iE£^ = k22lOD][H2] + k2o[OH][HD] - [Jja + Jn + (fei + k2s.)\H] + (^42 + k4i,)[He+] + fcep,4 + k,p,5)\HDO]. (5) 
at 

The photodissociation of HDO can either lead to OH + D ( Jto) or to OD -1- H (Jrt). The sum of the two photorates is the total HDO 
photodissociation rate and is simil ar to the H2O photodissociation rate Jg. The pho todissociation of HDO favors cleavage of the O-H bond 
over O-D bond with a 3 to 1 ratio Jshafer et al.lll989l ; IVander Wal et alj|l99dll99ll) . i.e. J-jh — SJia because the binding energy between O 



and D is stronger than that between O and H in water. In steady-state, the balance between formation and destruction leads to the ratio 

^ [HDO] ^ {k22[OD][H2] + k20[OH][HD])/{J7a + Jrt + (fc21 + k23)[H] + (fc42 + krs)[He+] + fcep.4 + fc,p,5) 

° [H2O] ks[OH][H2]/{J6 + ks[H] + k4i[He+] + kcp,3) ' 



_ [HDO] _ ( Dh,o \ ( k22([OD]/[OH])[H2] + k2o[HD\ \ 
° [H20\ \DnDo)\ h\H2] )' ^' 

where we define 

DH20 ^ J6 + k8[H] + k4ilHe+] + k,p,3 (8) 
and 

DhDO = Jva + Jib + (fc21 + k23)[H] + (^42 + k43)[He+] + fccp,4 + fc^p.S (9) 

to be the total destruction rates of H2O and HDO. The total UV photodissociation rates of H2O and HDO are similar Jzhang & Imr3l 19881) : 



Je ~ Jia + Jib- We further assume that rates with the deuterated species are similar to that of the main isotopologues: fcg = k2i + fc23, 
fc4i ~ ^42 + ^43, and kcp,4 + fccp,5 = kcp,3. If we can make the assumption that the total destruction rate of H2O and HDO are similar, 
i.e. DH2O — Dhdo, then the value of [HD0]/[H20] does not depend on the actual destruction mechanisms for H2O and HDO. The water 
fractionation /(HDO) is defined as 

f(Mnn\ [HD0]/[H20] ( Dh2o \ f k22 [OD]/[OH] , fc2(i^ /'.^^W'^22 feoX fc22 fc2o 

^ [HD]i[H2] = Vm^) [l^[HD]/[H2] + ^ j ^ Vm^) [^^^^^^ + " 1^^^^^^ + ^ 



(10) 

The ratio k2o / fcs is the ratio between the reaction of OH with HD to form HDO compared to the reaction of OH with H2 to form water. 
From the values in Table |4] the ^20/^^5 ratio is lower than I for gas temperature greater than 120 K. This term does not enhance the water 
deuteration fractionation. We can write analytically the ratio between rate fc22 and rate fcs : 

fc22 _ 1.55xl0-^^(r/300)^ «e-^««^/^ _ 756098fr/300)«°«e-^/^ (II) 
fcs " 2.05x10-12 (r/300)i-^2e-i660/T -U-75bU98(7/3UU) e . (11) 
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The ratio is weakly temperature-dependant. Since k-^^jk^ ~ 0.7—1 between 100 and 1000 K, water and hydroxyl radical deuterium fraction 
are similar at gas temperature greater than ~ 200 K: /(HDO)~ /(OD). Therefore water would be deuterium fractionated if OH is (i.e. 
/(0D)>1). 



2.2 Formation and destruction of OH and OD 

The water deuteration fraction /(HDD) is intimately linked to that OH /(OD). The rate limiting reaction for the formation of water is the 
formation of hydroxyl radical OH: O + H 2 — > OH + H with rate fci. T his reaction has a large energy barrier (_Ea=3163 K). The chemistry 
of OD (and OH) is described in details by Croswell & Dalgarnol 1 1985 ) for reactions without barrier. The production of OD mostly ( 



■ occurs 



through the rapid exchange reaction 

D + OH ^OD + H (12) 
once OH is present in the gas (reaction 16 with rate from Yung et al. 1988h . The forward reaction has no activation barrier but the reverse 



reaction (reaction 13) has a barrier of 810 K because OD is more stable than OH. Thus at T <1000 K, OD is favored. The role of OD for gas 
between 100 and 1000 K is similar to that of H2D^ for gas at T < 100 K for the deuterium enrichment. OD can also formed by the reaction 

HD + O^OD + H. (13) 

OD is destroyed by reaction with carbon ion 

C+ +OD ^CO + D+ (14) 
C"^ +OD ^ CO"^ + D, (15) 
and by photodissociation (with rate J5) 

OD + hu^O + D. (16) 
The steady-state abundance of OH and OD taking into account the most important formation and destruction reactions are 

. , ^ ki[0][H2] + {.h + ks.[H] + kii[He+] + kcp,z)[H20] + kiA[0][HD] + k22[OD][H2] + k-M[CO][H] + {Jja + k^p,5)[HDO] 

^ (fel5 + fcl6)[Z)]+fc4[i?]+fc5[H2]+fc44[He+]+fc6[0]+fc26[OD]+fc28[CO]+fc29[C]+fc3^ ^ ' 

,rjrr. k,[0][H2] 

^ ' k4H] + k5[H2] + k28[CO] + J4 + k,p,l + k29[C] + k39[C+] ^ ' 

After some algebra, we obtain the ratio 

[OH] ki[H2] 



(19) 



[O] J4 + fc4 [H] + k2S [CO] + kcp,l + k29 [C] + fc39 [C+] + fc44 [He+] + fccp.l ' 

which would be useful later in the analysis. For the deuterated hydroxyl radical at steady-state: 

. , ^ kj6[OH][D] + ki2[0][HD] + k2s[HDO][H] + k25[02][D] + k27[HDO][0] + k3-s[CO][D] + [Jjb + k,p,4)[HDO] 
^ ^ kn[H] + k22[H2] + k45[He+] + k24[O]+k26[OH] + k-34CO] + ka2[C] + k40[C+] + J5+k,p,2 

The rates are listed in Table 1 to 5. Neglecting the minor formation and destruction reactions (i.e. reactions with C and C"""), we simplify the 
ratio: 



[OD] ^ DoH / ki6[OH][D] + ki2[0][HD] + {Jrt + k2a[H] + k,pA)[HDO] \ 
[OH] ~ DoD I ki[0][H2] + (Je + k8[H] + kcp,a)[H20] 



(21) 



(22) 



where 

Dqh _ (fci5 + fci6)[D] + k4H] + ks[H2] + k44[He+] + fcefO] + k2ti[OD] + k2s[C0] + fc29[C] + fcagjC+l + Ji + fccp,i 

DoD ^ kir[H] + k22[H2] + k45[He+] + k24[0] + k2c[OH] + kai[CO] + hi2[C] + fc4o[C+] + J5 + fccp,2 

is the ratio between the sum of all destruction rate s of OH an d that of OD. Th e photodissociation rate of OH (rate J4) and OD (rate J5) are 
close I van Dishoeck & Dalgarndl 1984 : van Disho eck 1988; Croswell & Dalga rno 1985) (J4 ~ J5). From Table |4] reactions 5 and 22 have 



similar rates (fcs ~ ^22). Reactions 52 and 53 are ion-molecule reactions and we expect the rates to be of the same order of magnitude. Only 
reactions 4 (OH -1- H —*0 + H2) and 17 (OD + H ^OH -1- D) have very different rates, reaction 17 being much faster. But those reactions are 
important at low Av only where photodissociation dominates as destruction process. We assume that fee ~ fc24. 
The hydroxyl radical deuterium fractionation reads 



j(075) = [OD]/[OH] _ DoH ( k^(,(.[OH]/[0])[D] + k^2[HD] + {Jn + k23[H] + k 



cpa)Xd{{H20]/[0]) \ ( [HDl\ 

V[o]) J ' \ [H2] ) ■ ^''^ 



[HD]/[H2] DoD V ki[H2] + (Js + ks[H] + fccp,3)([H20]/[0]) 

The OD fractionation increases with larger amount of OH and atomic deuterium and decreases if oxygen is in the atomic form and 
deuterium is locked in HD. Large amount of OH allows the deuterium exchange reaction to occur. 

The knowledge of the abundances of H, H2, D, and HD is needed to estimate the hydroxyl deuterium fraction f{OD). 
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2.3 H and H2 

The formation of H2 occurs mostly on grain surface when the grain temperature is below 1000 K and the destruction is caused by photodis- 
sociation at low Ay and cosmic-ray/X-ray induced ionization and reaction with He^ in the UV free region. The temperature is low enough 
to avoid H2 destruction by atomic hydrogen. The steady-state balance between formation and destruction is: 

Ri{T„Ta)n„[H] = (/,,Ji + fec,4 + fc35[He+])[//2], (24) 

where fss is th e H2 self-shielding function dOraine & Bertoldill996h and Ri {Tg , Td) the molecular hydrogen formation rate on grain surfaces 
i TielenslboOSb 

i?i(r„rd) =4.4 X lQ-'^S{Tg,Ti) (25) 
where Tg and Td are the gas and dust grain temperature respectively and the sticking coefficient S is defined as 

S{Tg,Td) = ^ 10-2(rj + Td)i/2 + 2 X lO-STg + 8 X IQ-STl ' ^^^^ 

The sticking coefficient ensures that at high dust temperature atomic hydrogen does not stick onto grain surfaces. The number density of 
nuclei is 

na = [H] + 2[H2] ill) 
We obtain the atomic and molecular abundances: 

rrri ^ ( fssJl+kcA + k35[He+] \ ^ f fss.h+h,i \ 

' J \2R^nH + f.sJi + h,i + ki^[He+]) \2R^nH + fssJi + k^A J ^ ^ 

and 

r „ 1 f RiUh \ ^ / Rmn \ 

^ ^^^"^^ \2R^nH + fssJi+kQA + ks5[He+\) \2R^nH + fssJi+k^Aj ' 

Atomic and molecular hydrogen abundances are determined by the gas temperature, density, UV flux, extinction, the self-shielding function, 
and cosmic ray flux. 



2.4 D and HD 



While H2 is predominately formed on grain surfaces, formation o f HD can also occu r in the gas phase by reaction between D and H2 at 
high temperature and at density higher than about 5 x 10'' cm"'' 1 Le Petit et al.ll2002 ). HD is destroyed by photodissociation and reaction 



with atomic H while H2 is mainly photodissociated. H2 can self-shield against photodissociation, HD is shielded by dust only. Therefore 
deuterium remains atomic at higher extinction than H2. In hot and dense regions, an important formation route of HD is therefore: H2 + 
D — > HD -I- H (reaction 10 with a barrier _Ea=3820 K). Larger amount of atomic deuterium is needed to obtain high OD over OH ratio. 
The production of atomic deuterium is enhanced at high temperature via the conversion reaction HD -H H — > H2 -1- D. Another destruction 
mechanism of HD molecules involves photodissociation: HD + h.v ^\i + T>. Finally, reaction with H exchanges the atomic hydrogen with 
atomic deuterium: HD -1- H ^ H2 -1- D. At high Av, cosmic rays destroy HD. Reactions with H^ and atomic oxygen also destroy HD. The 
steady-state balance for [HD] then reads 

rrrr). ^ (fclQ [J/2] + ^2^1^ ) [jj] + fcsT [i/2] 

^ ^ J2+fcc,5+fcc,6+fcll[/i-]+fc36[i/+] + (fcl2+fcl4)[0]+fci8[0//]' 



where R2 is the formation rate on grain surfaces jLe Petit et a"i]|2002h density 



R2{Tg,Td) = &.i X io""s(r„rd) (^)'^' (31) 



and uh is the tot al number. This H2 form ation rate does not take into account grain chemisorption sites contrary to more sophisticated H2 
formation model 1 Cazaux & Tielenj|2002 ). We assume + fc^.s — ^C* ™d that the sticking coefficient is the same than for H2. 



We neglect the formation of HD via reaction with D^ and destruction via reactions with the protons, atomic oxygen and OH radicals. 
The steady-state balance leads to 

[D] ^ J2+fcc.4+fcll[iJ] 

[HD] kio[H2]+R2nH ^ ' 

If we can assume that most of the deuterium is locked in H and HD 

nD^[D] + [HD], (33) 
then we obtain HD and D abundances 
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\HD] = no ( ; ) (34) 



J2 + /cc,4 + A;ii[i?] + kio[H2\ + ii2nH 



and 



\D\=n ( J2 + fcc.4 + fcii[H] \ 

^ ^ ''\J2 + fcc,4 + fell [i/] + fcl„ [i/2l + i?2nij )' ^ ' 

At low Av, destruction of [HD] is dominated by photodissociation and at gas temperature a few 100 K, the formation of HD via reaction of 
atomic deuterium with molecular hydrogen is negligible 

[D] ^ _h_ 
[HD] RiUH' 

Even in obscured ( J2 ~ 0) and fully molecular regions {[H] ~ 0), cosmic-rays induced photodissociation and deuterium exchange reactions 
ensure that some atomic deuteriimi always remains in atomic form: 

[HD] - Rinn ^ ' 

The formation of HD on grains decreases dramatically for dust grain temperature above 100 K. 

2.5 Neutral and ionised Helium 

Where the steady-state abundance of ionised Helium is given by the expression: 

1 - fe5[//2]+fce-,3[e-] ^''^ 

Combining with the element conservation equation: 

UHe = [He] + [He+] (39) 
we obtain 

^ ^ k35[H2]+k,-,,[e-]-kc,3 ^ ' 

The He"*" abundance can be estimated only if the electron abundance is known. Our simple analytical analysis cannot provide an estimate of 
the He"*" abundance and we will omit in the rest of the paper reactions with He"*". 

2.6 Water fractionation ratio f{HDO) 

The results of the previous sections can be combined to derive an analytical formula of the [HD0]/[H20] ratio: 

■^o = ^^F^} = num/den, (41) 

[il2CJ 

where 

_ fc22 ( DoH Dh^o \ kui{[0H]/[0])[D]+kr2[HD] fcao [HD] 

"''"^ k^ \Dod Dhdo) ki[H2] + {JQ + ks[H]+kcp,z){[H20]l[0]) k5 [H2] ^ ' 

and 

ks \DoD DhdoJ ki[H2] + (Je + ks[H] + k,p-3){[H20]/[0]) ^ ' 

The abundances and abundance ratios in num and den have been derived earlier: 
[OH] ki[H2] 



[O] J4 + ki [H] + k2s [CO] + k29 [C] + fe39 [C+] + k44 [He+] + kcp,i ' 

[H2O] _ kr,[H2] 



(44) 
(45) 

[O] ~ (J4 + k4[H] + k28[CO] + fc29[C] + k39[C+] + k44[He+] + kcp,i){J6 + k8[H] + k4i[He+] + fccp.s) ' ^"^^^ 

The water fractionation is composed of 3 terms. At gas temperature greater than 200 K the second term ki2/ki dominates for both low and 
high extinction and the fractionation reaches 10-1000. Our analysis concerns high temperature gas only. In warm and dense region, the main 



[OH] Je + fcs [H] + k4i [He+] + k^p,s' 
and combining the two equations above, we obtain 

[H2O] _ klk5[H2] 
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Figure 1. Water deuteration fraction /(HDO). The upper panels show models with UV enhancement of factor 100 w.r.t. the standard interstellar' UV filed and 
for three increasing gas densities. The four curves correspond to the various dust extinctions. The lower panels are models with UV enhancement of 10*. 



oxygen carrier is water and that of carbon is methane CH4 and not CO at high Av and C at low Av- As state before, we further neglect 
reactions with He''' . 

At temperature below 100 K, the only fast chemical reactions are between ions and neutral species, which are not taken into account in 
our analysis. The starting point of water formation is H;^. Deuterium enrichment occurs via the deuterated equivalent of Hjj', H2D'''. H2D+ 
has a lower zero-point energy, which favors the deuteration reaction H|j' + HD ^H2D+ + H2. The rate of the back reaction has an energy 
barrier of 230 K. 

The water deuteration fraction enhancement is caused by the fast deuterium exchange reaction with OH (D+OH ^OD + H). The inverse 
reaction has an energy barrier of 717 K. Significant amount of atomic deuterium is possible at low and intermediate extinction. We plotted 
in Fig. [T] the ratio [HD0]/[H20] as function of the gas temperature for an impinging UV enhanced by a factor 10* and 10'^, for different 
extinctions (Av=l, 5, 10, 15), and for two gas densities (nn = 10* and 10^^ cm~"^). Disc surfaces around young accreting T Tauri stars 
receive ~ 10* to lO'' times the amount of standard interstellar UV. The UV field, extinction, density, and gas temperature values are typically 
found in regions of discs where water is abundant. In the following section, we will use a time-dependent thermo-photochemical code to 
support our assumed values for the parameters. The figure shows that large /(HD0)=[HD0]/[H20] ratios (up to a few 100) are possible 
for gas temperatures up to 500-600 K and up to extinction of 10, although the average value for / is lower than 10. At T >500, the back 
reaction OD-l-H — >OH-l-D starts to destroy efficiently OD. The figure also shows the large range of values for /(HDO) (0.1-10^). The /(HDO) 
curves testify of the sensitivity of the abundances on the gas temperature, which reflects the exponential nature of the Arrhenius' law for 
neutral-neutral reactions. 
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Table 1. Stellar and disc parameters 



Stplliir inn';'; 




S M^ 


Stellar luminosity 


r 


n 7 T ^ 


Effective temperature 


T a 
erf 


4400 K 


Disc mass 


Md 


10"^ M© 


Disc inner radius 




0.1 AU 


Disc outer radius 




300 AU 


Vertical column density power law index 


e 


1 


dust to gas mass ratio 




0.01 


dust grain mateiial mass density 


Pdust 


2.5gcm-3 


minimum dust particle size 




0.01 lira 


maximum dust particle size 




100 /im 


dust size distribution power law 


P 


3.5 


Cosmic ray flux 


CR 


1.7 X 10-17 s-l 


ISM UV filed w.rt. Draine field 


X 


1.0 


abundance of PAHs relative to ISM 


/pah 


0.1 


a viscosity parameter 


a 


0.0 



3 TIME-DEPENDENT MODELLING 



3.1 ProDiMo photochemical code 

We modelled the chemical abundances of a 10~^ Mq disc around a typical TTauri star (Teff 4400 K, log g=A.Q, Z=1.0) using the pho- 
tochemical code ProDiMo . ProDiMo combines frequency-dependent 2D dust-continuum radiative transfer, kinetic gas-phase and UV 
photochemistry, ice formation, and detailed non-LTE heating and cooling balance. The major improvement over previous studies is that the 
density structure is determined by the gas pressure that is computed by detailed chemistry and energy balance d and not by assum i ng tha t 
the gas and dust have the same temperature. Detailed description of the code are given by Woitke et al. ( 2009 ) and in Kamp et al. ( 20091) . 
The code has been used to determine the water abundance in the inner disc around a typical Herbig Ae star "l lWoitke et alj|2009l) . The most 
recent additions to the code includes PAH chemistry, time-dependent chemistry, deuterium chemistry, and generation of Spectral Energy 
Distribution and spectral lines. The last two features are not used in the modellin g performed for this paper. Table [T] summarized the input 
parameters to the mode l. The stellar spectrum w as generated using PHOENIX (Bro tt & Hauschild3 l20051 with the addition of chromospheric 
flux from HD 129333 jDorren & Guinanlll994f) . Although the model extends to 300 AU, we focus here only on the inner 3 AU as we are 
interested in the [HD0]/[H20] ratio in the terrestrial planet forming region of discs. 

The chemical netw ork includes a total o f 187 deuterated and non-deuterated gas and ice species. Most reaction rates are taken from 
the UMIST DATABASE JWoodall et alj|2007h . A dditional reaction rat e s were compiled from th e NIST chemical kinet ic da tabase. The rates 
involv ing deuterated species are described in e.g., Roberts et al. 1 2004 ), Roberts & Millai 1 200d) . Charnlev et al. 1 1997 ), and Brown & Millai 
1 19891) . Species can freeze-out onto grain surfaces and desorb thermally or upon absorption o f a cosmic-ray or a UV ph oton (photodesorption). 
Grain surface reactions were omitted apart from the grain surface formation of H2 and HD 1 Cazaux & Tielensll2002 ). The photodissociation 
cross-sections are taken from the LEIDEN DATABASE described in lvan Dishoeck et al. 1 200a) ^ 

The chemical abundances in the disc were established in three stages. First ProDiMo determined the UV field, hydrostatic density, 
gas and dust temperature, and chemical abundance structure self-consistently assuming steady-state chemistry. Second we computed the 
chemical abundances of gas and solid species for a 1 Myr old molecular cloud with density of 5 x 10* cm~^ and gas and dust temperature of 
15 K from diffuse cloud initial abundances, where the elements are in neutral or ionized atomic form (see Table|2j. Third we ran ProDiMo 
in the time-dependent chemical mode to simulate the chemical structure of a 1 Myr old disc using the results of the molecular cloud run as 
initial chemical abundances and the disc properties computed in the first stage. All other disc properties (UV field, density and temperature 
structure) were fixed in this last stage. 

The three-stage method mimics the incorporation of molecular cloud materials and their subsequent evolution in the disc. It also makes 
it possible to comp are the chemical abundances obtained with the time-dependent model and at steady-state. Our approach differs from 
who solve the chemistry in a Lagrangian frame and their disc evolves according to viscous spreading. However their 



Visser et al. 



3 comp a 
1 2OO9I) \ 



number of species and reactions are limited. 



3.2 Model results 



Figure [2] show the UV field, density, and gas and dust temperature structure in the inner 3 AU. The disc structure is discussed in lWoitke et al 



1 20091) . In this paper we focus on the water and HDO abundances. Figure|3]show the water and HDO abundances at steady-state on the left 
and for a disc of 1 Myr old on the right. Apart from the H2O HDO in the midplane beyond 1.5-2 AU, the water and HDO abundances 
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Table 2. Typical diffuse cloud abundances used as initial abundances for the molecular cloud chemical calculation. Species with ionization potential (IP) 
higher than 13.6 eV are neutral while species with IP below 13.6 eV are ionized. The other species have negligible initial abundance. 



Species 


log(n(X)/nH) 


H 


0.0 


D 


-5.0 


He 


-1.125 


C+ 


-3.886 





-3.538 


N 


-4.67 


S+ 


-5.721 


Si+ 


-5.1 


Mg+ 


-5.377 


Fe+ 


-5.367 


PAH 


-6.52 



Table 3. Mass and average gas and dust temperature for the three gas-phase water and HDO locations. The paper focuses mostly on the inner midplane. 



Location Mass < Tgas > < Tdust > 

(Mo) (K) (K) 



H2O 



Inner midplane 


2.6 


X 


IQ- 


-8 


178 


177 


Cold belt 


1.1 


X 


10" 


-9 


16 


16 


Hot layer 


2.1 


X 


10" 


-12 


1147 


120 



HDO 



Inner midplane 1.2x10-10 140 140 
Cold belt 3.8 X 10"" 16 16 

Hot layer 2.6x10-1" 474 102 



are similar for the steady-state and time-dependent chemistry models. The differences stem from the fact that at steady-state even very slow 
reactions will impact the chemistry. In this case, the dust temperature is low enough (Td < 150 K) such that all water and HDO should be 
frozen-out. In the time-dependent model, significant amount of water and HDO remains in the gas-phase in the inner 1 AU or above the ice 
zone. The differences between the steady-state and time-dependent abundances occur in low gas-phase water abundance regions and thus do 
not impede our analysis of deuterium enrichment for gas-phase water either outside the water freeze-out zone. 

The [HD0]/[H20] ratio for the time-dependent model is plotted in Figure|4] The water abundance levels at 10~® and 10~** are overlaid. 
Water is abundant as soon as hydrogen is in molecular form. In our model, water and HDO are frozen onto grain surfaces in the midplane 
beyond ~1.5 AU as shown in Figure. [5] The gas-phase [HDO]/[H20] ratio decreases with radius. At 0.5-1.5 AU in the mid-pl ane, the ratio 



is a fe w 10 ^-10 ^ . Gas-phase water and HDO are located in three zones: the inner midplane, the cold belt, and the hot layer jWoitke et al 
2009h . The mass and average gas and dust temperature are given in TablejS] 



Most of the gas-phase water (n(H2O)/ nH=10~^-10~^) is located in the inner plane close to the star {R <l AU). In the inner 30 AU a small 
fraction is found above the mid-plane 1 Woitke et alj|2009|). where the vertical extinction Ay is between 1 and 10. Although of insignificant 
mass, the molecules are hot and emit strongly (e.g. Carr & Naiital 20081 : Salvk et al. 200 sh . In the outer disc, gas-phase water is found in 
a col d belt, which is san dwiched between Ay ^ I and 5 where efficient photodesorption maintains some water molecules in the gas phase 
(see jWoitke et al]|2009h ). Since HDO and H2O have similar adsorption energy, the gas-phase abundance of both species are co-located (see 
FigO. We plotted the vertical column densities for H, D, OH, OD, H2O, HDO, H20#, and HDO# in Fig.|6] This fi sure shows that the 
column density ratios HDO/H2O, HD0#/H20#, and OD/OH stay relatively constant in the inner disc, with values much higher than the 
elemental D/H ratio of 10^'^. 

In the high water abundance region, the [HD0]/[H20] ratio is higher than lO^'^. Using the mass of H2O and HDO in the inner midplane 
listed in Table[3] we derived an average [HD0]/[H20] ratio of 4.6 x 10^'^, which is 30 times higher than the Earth Mean Ocean Water value 
(~ 1.49x 10~^) and close to cometary values. The actual [HD0]/[H20] values in our models depend on the disc parameters. Future studies 
with focus on the effect of disc properties (disc mass, radius, ...) on the [HD0]/[H20] ratios. 

Overall, the [HD0]/[H20] ratios in the inner disc midplane by the time-dependent code ProDiMo are consistent with the analytical 
results. Our results contradict the simple decrease in [H20]/[HD0] with increasing temperature if thermochemical equilibrium ratios are 
assumed. 
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Figure 2. The four panels show the gas density n < h > (upper left), strength of the UV field with respect to the Draine field x (upper light), the dust temperature 
T(j (lower left), and the gas temperature Tg with dust and gas temperature overlaid (lower right) for the inner 3 AU computed by ProDiMo. The structure is 
consistent with the chemical abundances. 



4 DISCUSSION AND CONCLUSION 

A simple analytical analysis shows that high water fraction is possible in a simple deuteration chemical network for gas hotter than 100 
Kelvin if n eutral-neutral reac t ions a re included. The analytical results are supported by the time-dependent photochemical model results. In 
their studv lWillacv & Woodsl ( 120091) also found that HDO can be abundant in the inner disc. 

The water deuterium fractionation is determined by the ratio between the rate of formation of OD via O + HD and the rate of for- 
mation of OH via O + H2. Lower zero-point energy for HD makes the first reaction faster than the second one above 200 K. At low 
temperature (T < 100 K), fast deuteration of H;^, the main driver of cold gas-phase chemistry ensures that main molec ules are deuterium- 
enriched. Another reason for higher HDO/H2O is the preferential branching into OD + H when HDO is photodissociated dShafer et al .1 19891 : 
Vander Wal et al.lll99dll99ll) . 

From our work, it is also theoretically possible to have high water deuterium fractionation from gas-phase photo-chemistry above 200 K. 
Deuterium errrichment can also occur at high temperature because of the energy difference in activation barrier for deuterium exchange and 
the back reactions. Finally, water deuterium enrichment is not as stringent a constraint as thought for the origin of water on Earth. Deuterium- 
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Figure 3. Gas-phase H2O and HDO abundance in the inner 3 AU computed by ProDiMo. The left panels show the steady-state abundances while the right 
panels show the time-dependent abundances for a 1 Myr old disc. 



enriched water may have been synthesised at '--^ 1 AU and incorporated directly onto sihcate dust grains jStimpfl et allbood) . Water could 
have been incorporated on Earth already at the stage of planetesimals accretion. However our simple analytical steady-state chemical model 
and the 2D modelling with ProDiMo do not include effects of turbulent mixing in protoplanetary discs, both radial and vertical. 
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Table 4. Principal reactions 





Reaction 




A 

(cm'^ 


s" 






B 


(K) 


Reference 




O + H2 


OH + H 


3.14 


X 


10" 


■13 


2.7 


3150 


UMIST 




OH + H2 - 


-* H2O+ H 


2.05 


X 


10- 


12 


1.52 


1660 


UMIST 


fcio 


H2+D-* 


HD + H 


7.50 


X 


10- 


■ 11 




3820 


Zhang & Millar 1989 


fell 


HD + H- 


*H2 +D 


7.50 


X 


10- 


■ 11 




4240 


Zhang & Millar 1989 


ki2 


+ HD- 


»OD + H 


1.57 


X 


10- 


■12 


1.7 


4639 


Joseph, Trahlar & Ganet 1988 


kie 


OH + D- 


»OD + H 


9.07 


X 


10- 


-11 


-0.63 




Yung et al. 1988 


k20 


OH + HD 


HDO + H 


0.60 


X 


10- 


■13 


1.9 


1258 


Talukdar et al. 1996 




OD + H2 - 


HDO + H 


1.55 


X 


10- 


■12 


1.6 


1663 


Talukdai- et al. 1996 




Reaction 




q 








7 




Reference 



(s-i) 



J2 HD + hu^U + D 2.6x10"" 2.5 - Le Petit et al. 2002 

J4 OH + hi/^O + H 3.5x10-"' 1.7 - van Dishoeck 1988 

Ref. UMIST: IWoodall et alJiioOTi) 



Table 5. Neutral-neutral and radical-neutral reactions 





Reaction 




A 








B 




Reference 








(cm'^ 


s" 






(K) 


(K) 




Ri 


H-hH-i-M 


^H2 










- 


- 


see text 


H.2 


H-hD-i-M 


-»HD 














see text 


ki 


0-HH2 ^ 


OH-hH 


3.14 


X 


10- 


-13 


2.70 


3150 


UMIST 


k2 


O2 -hH^ 


OH-hH 


2.61 


X 


10- 


10 




8156 


UMIST 


fcs 


O2 -1- H2 - 


+ OH -H OH 


3.16 


X 


10- 


■10 




21890 


UMIST 


^4 


OH-hH^ 


0-HH2 


7.00 


X 


10- 


-14 


2.80 


1950 


UMIST 


k5 


OH + H2 - 


-* H2O+ H 


2.05 


X 


10- 


12 


1.52 


1660 


UMIST 


ke 


OH -1-0^ 


O2 -f H 


1.77 


X 


10- 


■11 




-178 


UMIST 


k7 


OH -H OH - 


-^0-l-H20 


3.87 


X 


10- 


■13 


1.69 


-469 


UMIST 


ks 


H2O-1-H- 


^ OH -1- H2 


1.59 


X 


10- 


■11 


1.20 


9610 


UMIST 


kg 


H2 0+0^0H-(-0H 


1.85 


X 


10- 


■11 


0.95 


8571 


UMIST 


kio 


H2-l-D^ 


HD-f H 


7.50 


X 


10- 


■11 




3820 


Zhang & Millar 1989 


kii 


HD-hH^ 


H2 -l-D 


7.50 


X 


10- 


■11 




4240 


Zhang & Millar 1989 


ki2 


0-l-HD^ 


OD-l-H 


1.57 


X 


10- 


■12 


1.7 


4639 


Joseph, Truhlar & Garret 1988 


ki3 


OD-hH^ 


0-1- HD 














ki3 = ^4 is assumed 


fcl4 


0-i-HD^ 


OH -l-D 


9.01 


X 


10- 


13 


1.9 


3730 


Joseph, Truhlar & Garret 1988 


ki5 


OH-hD^ 


O + HD 














ki5 = k4 is assumed 


ki6 


OH-l-D ^ 


OD-t-H 


9.07 


X 


10- 


■11 


-0.63 




Yung et al. 1988 


ki7 


OD + H^ 


OH + D 


1.26 


X 


10- 


■10 


-0.63 


717 


Yung et al. 1988 


kis 


OH -H HD - 


-^H20-l-D 


2.12 


X 


10- 


13 


2.7 


1258 


Talukdar et al. 1996 


ki9 


H2O-1-D - 


OH -1- HD 














fcig = fcg is assumed 


k20 


OH + HD - 


^ HDO + H 


0.60 


X 


10- 


■13 


1.9 


1258 


Talukdar et al. 1996 


k2i 


HDO + H ^ OH -1- HD 














fc2i = 0.5 X kg is assumed 


k22 


OD + H2 - 


HDO + H 


1.55 


X 


10- 


12 


1.6 


1663 


Talukdar et al. 1996 


k23 


HDO -1- H - 


OD -1- H2 














k23 = 0.5 X kg is assumed 


k24 


OD-fO^ 


O2 -hD 














k24 = ks is assumed 


k25 


02-hD^ 


OD-hO 














k25 = k2 is assumed 


k26 


OD -1- OH - 


-H HDO 














k26 = kr is assumed 


k27 


HDO + 


^ OD -t- OH 














k27 = fcg is assumed 


k2S 


OH + CO - 


CO2 + H 


1.17 X 10- 


13 


0.95 


-74.0 


UMIST 


k29 


OH-l-C^ 


CO + H 


1.10 X 10- 


■10 


0.5 


0.0 


UMIST 


kao 


CO-^H^ 


OH-l-C 


1.10 X 10- 


■10 


0.5 


77700.0 


UMIST 


kai 


OD + CO - 


CO2 + D 














ksi = k28 is assumed 


k32 


OD-l-C-> 


CO-hD 














^32 = k29 is assumed 


k33 


CO-hD^ 


OD-hC 














^^33 = ^^30 is assumed 



Ref. UMIST: IWoodall et alJiioOTi) 
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r [AU] 

Figure 4. [HDO]/[H20] in the 3 AU of a 10~^ disc as computed by the photochemical code ProDiMo . The contours indicate the regions where gas-phase 
water abundance is 10~® and 10~*. 
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Figure 5. H2O and HDO ice abundance in the inner 3 AU. The dash-line indicate the region where the dust temperature is 150 K. H2O and HDO freeze at 
~150 K. 



Table 6. Photodissociation and photoionisation reactions 



Reaction q 7 Reference 





Jl 


U2 +hi/ — 


»H + H 


3.4 X 10" 


ll 


2.5 


self-shielding factor, see text 


J2 


HD + hu- 


^ H + D 


2.6 X 10- 


11 


2.5 


LePetit et al. 2002 


J3 


CO + hu- 


^C + 


2.0 X 10- 


10 


2.5 


UMIST 


Ja 


OH + hu- 


-* O + H 


3.5 X 10" 


10 


1.7 


van Dishoeck 1988 


Js 


OD + hu- 


^ O + D 


4.0 X 10" 


10 


1.7 


Croswell & Dalgamo 1985 


Je 


H2O + hu 


^H + OH 


5.9 X lo- 


10 


1.7 


UMIST 


Jra 


HDO + hu 


^OH + D 








Jra = 0.25J6, see text 


Jib 


HDO + hu 


^OD + H 








Jrb = 0.75J6, see text 


J9 


C + hu^ 


C+ -l-e 








UMIST 



Ref. UMIST: IWoodall et alJiioOTi) 



Table 7. Ion-neutral reactions 





Reaction 


k 

(cm^ s~^ 


) 


Reference 




H+ -H D ^ D+ + H 


1.0 X 10 


-9 


Watson (1976). UMIST 


^35 


He+ -1- H2 ^ H + H+ + He 


1.1 X 10 


-i3(T/300)-0-24 


UMIST 


^36 


HD -1- H+ ^ D+ -H H2 


1.0 X 10 


-9 g-464/T 


UMIST 


^37 


H2 -1- D+ ^ H+ -l-HD 


2.1 X 10- 


-9 


UMIST 


fc38 


H -1- D+ H+ -H D 


1.0 X 10 


-9 


Watson (1976). UMIST 


^39 


OH -H C+ -+ C0+ -1- H 


7.7 xlO- 


10 


Prasad & Huntress (1980), UMIST 


fc40 


OD -H C+ ^ CO+ -(- D 






k4o = ksg is assmued 


fc41 


He+ -1- H2O ^ OH -H He -H H+ 


6.0 X 10 


-11 


UMIST 


^42 


He+ -1- HDO ^ OD -1- He -1- H+ 






A;42 = 0.75fc4i is assumed 


fc43 


He+ -1- HDO ^ OH -1- He -1- D+ 






^43 = 0.25fc4i is assumed 


^44 


He+ -1- OH 0+ -H He -H H 


1.1 X 10 


-9 


UMIST 


fc45 


He+ -1- OD ^ 0+ -H He -H D 






/C45 = fc44 is assmued 



Ref. UMIST: IWoodall et alJiioOTi) 
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Figure 6. Half disc vertical column density for vaiious species in the inner 3 AU. 
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Table 8. Cosmic-ray induced ionisation and recombination reactions = 5 X 10 s ^). 





Reaction 


k 




Reference 




T-T J- rr * 4- p~ 


0.46 X ( 




UMIST 




D + cr — » D+ + e~ 


0.46 X C 




UMIST 


kc,3 


He + cr ^ He+ + e~ 


0.5 X C 




UMIST 


kc,4 


H2 + cr ^ H+ + H + e" 


0.04 X C 




UMIST 




HD + cr ^ H+ + D + e" 






fc^ 5 = 0.5fc^_4 is assumed 




HD + cr ^ D+ + H + e~ 






fe^ 5 = 0.5fc^_4 is assumed 


K~,l 


U+ +e~ ^H + hu 


3.5 X 10-i2(T/300)-0-^o 




Prasad & Huntress (1980) 


^e- ,2 


D+ +£- ^D + hu 






fcg- 2 = k^- ^ is assumed 


^e- ,3 


He+ + e~ — > He + hi/ 


4.5 X 10-i2(T/300)-'' «^ 




Prasad & Huntress (1980) 


^cp, 1 


OH + CRPliot ^ + H 


1.3 X 10-i''{509/(l - w)) 




UMIST 


^cp,2 


OD + CRPliot ^ O + D 


1.3 X 10-1''{509/(1 - w)) 




kcp,2 = fccp.i is assumed 




H2O + CRPhot ^ OH + H 


1.3 X 10-i''{971/(l - w)) 




UMIST 




HDO + CRPliot ^ OD + H 


0.75 X 1.3 X 10-1^(971/(1 


-w)) 


UMIST 


^cp,^ 


HDO + CRPliot ^ OH + D 


0.25 X 1.3 X 10-1^(971/(1 


-w)) 


UMIST 
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